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I present a theory of topological heavy-fermion semiconductors based on the large-A' symplectic representa- 
tion for the electron spin. The theory is exact in the limit when the number of spin flavors N = 2kis infinite. 
I find that both weak and strong topological insulating phases exist for k < 3. Furthermore, for k > 3 the 
weak topological insulating state fully suppressed while only strong topological and trivial insulator states sur- 
vive. In addition, using the mean-field theory results, I consider the tunneling into topologically trivial and 
non-trivial phases of a generic heavy-fermion insulators by calculating the differential tunneling conductance. 
The implications of the presented results for the existing heavy-fermion semiconductors are discussed. 
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I. INTRODUCTION 

Heavy fermions are a class of complex materials in which 
the effective mass of conduction electrons greatly exceeds 
the bare electron mass due to the strong hybridization be- 
tween the conduction and localized /-electronJ02. The prop- 
erties of heavy-fermion materials present an interest for re- 
searchers from the perspective of both technology and fun- 
damental science. While the technology-driven research on 
heavy fermions is mostly focused on problems related to 
energy conservation and storage, fundamental physics re- 
search takes two interconnected directions: the investigation 
of the microscopic origins of unconventional superconductiv- 
ity (with CeCu2Si2, UBei3 and UPtj as first examples of this 
phenomenon in solid state systerr!"^, and discovery and anal- 
ysis of the novel states of quantum matter. 

Despite many years of experimental and theoretical re- 
search, the microscopic mechanisms responsible for the emer- 
gence of various quantum states in these materials remain 
unclear The hidden-order phase in URu2SiJS^, rnagnetic 
field-induced non-Fermi-liquid behavior in YbRh2Si?H! low- 
temperature metallic conductivity in heavy-fermion semicon- 
ductor SmBe and Ce-based compounds'^ superconducting 
response to Yb -doping in CeCoInj'^ La- and Y-doping in 
CeCu2Sij22llIl and recently discovered quantum critic ality in 
y6-YbAlB4^ are just a few examples. All these states emerge 
via a physical process unique to heavy-fermion materials in 
which strong interactions between conduction and localized 
/-electrons operate in an environment of very strong spin- 
orbit coupling. 

SmBs is a pro totypical example of a heavy-fermion 
semiconducto^^^'. Interaction between Sm p- and tZ-orbital 
conduction electrons and localized /-electron states leads to 
an opening of the hybridization gap at T* ^ lOOK. The av- 
erage electron /-level occupancy is well below one, rtj < 
1, demonstrating the strongly mixed valent nature of this 
material^^ Transport measurements in SmBe show an in- 
crease in resistivity below T ^ 50K and then saturates at very 
low (below 5K) temperatureJ^i'f^ll xhe value of residual re- 
sistivity grows when the quality of the sample increases and 
is of the order of psa, ~ 30 Q-cm. This value is incompatible 
with the one originating from the metallic conduction in the 
presence of disorder induced scattering. The origin of low- 



temperature conductivity still re mains poorly understoo d'^ l ^^ l 
however, it was recently proposecP^'^ll that SmB5 is a topolog- 
ical insulator and finite low-temperature conductivity can be 
due to topologically protected metallic surface states at sam- 
ple boundaries. 

There are two major theoretical challenges in understand- 
ing the anomalous transport properties of SmBf,. The first 
challenge is due to the strong Hubbard interaction between 
the conduction and /-electrons. The second challenge is that 
the symmetry of the lowest lying crystalline field multiplets, 
which are hybridized with conduction electrons and determine 
the symmetry of the hybridization gap, is not known. These 
challenges make the formulation of full microscopic transport 
theory in SmBe quite challending. Nevertheless, some general 
features of the heavy-fermion semiconductors in regards to the 
topological features in their band structure can be described 
on a more general (i.e. less material dependent) level. One 
important question, for example, is the question of stability of 
weak and strong topological insulating phases depending on 
degeneracy of the local /-level. 

The basic model which is thought to capture the main as- 
pects of the physics of the heavy-fermion semiconductors is 
the Anderson lattice model. The Hamiltonian can be written 
as a sum of the following three terms: 

H^H,+Hf + Hh. (1) 

Here the first term He describes the conduction electrons 

Hc = Yj ^k4o-'^k<r, = - ^ ^ cos ki - fi, (2) 

k,tr i=x,y,z 

where cr denotes electron's spin projection, t is the hopping 
amplitude (equal to bandwidth) and fj.^ is the chemical poten- 
tial. Consequently, /-electrons are described by 

j a=±\ ia (ij) a=±l 

where fj^ creates an /-electron on site j in a state a of a lowest 
lying multiplet A^p-degenerate multiplet denoted by F (below 
we consider Kramers doublet only, so Nr - 2), e^°' is the /- 
electron energy and U > Ois the Hubbard interaction between 
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the /-electrons. The last term in (j3]l yields a very weak hole- 
like dispersion for /-electrons to enforce the fully gapped in- 
sulating state. We emphasize that index a is not a spin index 
due to the presence of the strong spin-orbit coupling. Lastly, 
the third term in Eq. ([T]) accounts for the interaction between 
conduction and the / electrons 



j,a=±\ 



(4) 



In Refs. II33I34I the phenomenological analysis of the An- 
derson lattice model based on the low-energy expansion of 
the /-electron self-energy has been used to analyze the 
topological structure of the resulting heavy-fermion semicon- 
ducting state. In this paper, I will resort to a microscopic ap- 
proach based on the large-A^ slave-boson theory and analyze 
the topological structure of the insulating state. I employ the 
symplectic SP(A^) {N - 2k,k - 1,2, ...) representation for the 
electronic operators to properly describe time-reversal sym- 
metry ofthe electronic states. In agreement with the previous 
resultPl211 I find that for = 2 and N - A there appears two 
(weak and strong) topologically non-trivial states depending 
on the relative position between the renormalized /-level and 
the chemical potential of the conduction band. Moreover, I 
found that for the large value of > 4 there is only strong 
topological insulating state. In addition, I will discuss the 
tunneling into a topologically non-trivial heavy-fermion semi- 
conductors. 

In the next Section I will present the large-A^ mean-field 
theory of topological heavy-fermion semiconductors, and the 
calculation of the tunneling conductance of the surface states 
in weak topological insulator. The Section II is devoted to the 
discussion of the results and conclusions. 



II. SYMPLECTIC SLAVE-BOSON THEORY 

Large-A^ slave-boson mean-field theories (A^ is the degen- 
eracy of the /-electron level) utilize the naturally small pa- 
rameter 1/A^ to determine the thermodynamic properties of 
heavy-fermion materials by expanding near exactly solvable 
limit^'^^"* of A^ — > oo. Recently, a novel large-A^ expansion 
methods have been developed to account for the specific sym- 
metries of the problem (see |45 1 and references therein). In 
what follows, we generalize our model from SU(2) symme- 
try group to SP(A^) with A^ = 2k, so that the spin summation 
run over k spin indices, a,cr e {+\,+k]. The importance of 
using the SP(A^) subset of SU(A^) group clearly lies in the re- 
quirement for the proper description of the states related by 
time-reversaF^. In the SP(A^) version of the theory, the form- 
factor matrix acquires a block-diagonal form of identical 2x2 
blocks. For the subsequent saddle-point analysis we find it 
more convenient to use the path integral formulation. 

The slave-boson approximation corresponds to (i) taking 
the limit t/ — > oo, which corresponds to projecting out the 
doubly occupied states and (ii) introducing the constraint 



which guarantees the local moment at the /-site, i.e. nj - \: 



(5) 



The partition function corresponding to the model Hamil- 
tonian ([TJ with constraint condition (|5| reads : 

n//3 



J — J D(b,b\f,f\c,c^)exp - J L(T)dT 

-tt/B V 



where Lagrangian L{t) is 



, (6) 



;' i I ty= 1 L \ / / 



k a=\ 



+ Z Z ([<l'rk].,.'^«'/>,ck, +h.c.) 

'^■'^ i,k a.p=l 

f N \ 



d 

d^ 



\a=l 



(7) 



Here we use the same notation for the form-factor matrix, al- 
though now it has a block diagonal form of k(= N/2) blocks 
each of dimension 2x2. We have also introduced the field 
Aj to enforce a constraint. Finally, we have rescaled the hy- 
bridization amplitude V V/ y/N for the bookkeeping pur- 
poses (see below). 

Now we can integrate the conduction electrons by making 
the following transformation 



cl — > cl 



V 



Z £[4'rk]<./se-'''«'(5, + ^kr'hb].^ 



i 0=1 

N 



(8) 



V 



2 Z['^rkWe*'*'/>,(5r + ^k)"' 



In what follows, it is convenient to write the 2x2 form-factor 
matrix as follows 

Ork = Mnk -A 0k = ^Tr[<l)[^<l)rk] (9) 
where t are Pauli matrices and Wk is a unit vector. We obtain 



I I J a= 1 ^ 

k a=l ^ ' > \a=\ ) 



N 



2 £ A„,(k)."'*-«^)/.>,(5.+^k)-'/;,/^], 



(10) 
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where we have introduced 



A. mean-fleld solution 



1 ^ 



The action with the Lagrangian above is quadratic in 
fermionic operators, which can be integrated out to give an 
effective action in terms of the slave fields only. Since the 
resulting expression is quite cumbersome we will not give it 
here. Instead, we proceed with the saddle-point analysis. 
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FIG. 1: Phase diagram found from the solution of the slave-boson 
mean-field equations l[T6j. Kondo hquid state corresponds to the sit- 
uation when the slave-boson amphtude a = 0. I have used the fol- 
lowing values for the input parameters: = O.lf, = -1.05/. 
When the number of the fermionic flavors exceeds four. A' > 4, there 
exists only strong topological insulating phase and weak topological 
insulating state disappears. 



Mean-field (saddle-point) approximation corresponds to 
the following values of the bosonic fields: 

^(t) = VA^a5,,o, i/lq(T) = (e/ - sf)6,,o, ( 1 1 ) 

where both a and s/ are r-independent. Now, we can use the 
Matsubara frequency representation and integrate out /-fields, 
since the action is quadratic in these fields. These yields: 



-nIP 



Seff^N{ef-ef)(a^-qN) 

- 2NT y y log[(/w - coikXioJ - co2k)]. 



(12) 



/oj k 



where we have introduce the parameter - jj- Moreover, 
functions wi_2k describe newly formed energy bands 



Wl,2k 



^k + £/k ± 7(^k-£/k)2 + 4(yfl<^k)2 

Ef^-Ef-\-Nah\^, h\^- -ty ^ cos/:/. 



(13) 



We note that the newly formed band spectrum corresponds to 
the eff'ective Hamiltonian 

The reason we invoke the effective Hamiltonian is that it will 
allow us to analyze the topological structure of an insulating 



state. In Eq. ( 14 1 ]_ denotes the unit 2x2 matrix. To determine 
the parameters a and ej we, of course, have to minimize the 
effective action. In addition, we have to keep in mind that the 
total number of electrons must be conserved. Specifically, for 
an insulator, we have to require that we will have one conduc- 
tion electron per /-electron, so that 

ric+rif-N. (15) 

Minimization of the effective action together with the con- 



dition for an insulator ( 15 i yields the following system of the 
mean-field equations: 

{Sf- s'p)a + tYj [Nahi,Aff(k, icj) + V(pkA /^(k, ioj)] = 0, 

/t(>,k 

(a^-qj,) + Tj^Yj^ ff(k, icj) = 0, 

ioj k 

{qN - a^) + r ^ ^ A„.(k, iu) = 1, 



iw k 



(16) 
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where the functions Aaft(k, ito) are defined by 
Affikjcj) = 
A /c(k, io)) = 













yfl0k 






- ^k)('w - -E/k) - 






/ W - ^k 






- Efk)(icj - Efi,) 





(17) 



l'^^, where 6 define the di- 



To solve ( [T6] l we still need to specify the momentum de- 
pendence of the hybridization gap, (p^. In what follows, we 
adopt the choice of the form-factors from Refs. |33 34| and 
consider function 0k which at small momenta k is <p^ - 

^ ^[12cos(20) + 5(3 -1- cos(46'))]i 

rection of the unit vector k, associated with the point on the 
Fermi surface. 

To analyze the topology of the bands governed by the ef- 
fective Hamiltonian (14 1 we use the fact that topology is in- 
variant under any adiabatic deformation of the Hamiltonian. 
We begin our study with a tight-binding model for a Kl on 
a simple cubic lattice. Our choice of hybridization ensures 
that the mean-field Hamiltonian (Eq. 14 1 is a periodic func- 



tion satisfying '?4//(k) - 'K(.//(k + G). The technical analysis 
is readily generalized to more complicated cases as discussed 
below. The most important element of the analysis is the odd 
parity form factor of the / electrons, <l'r(k) = -(l)r(-k). This 
parity property together with the absence of the nodes in the 
hybridization gap across the Brillouin zone (BZ) are the only 
essential input as far as the topological structure is concerned. 

To evaluate the topological indices we use the results of 
Fu and Kane [46] who have demonstrated that in an insu- 
lator with time-reversal and space-inversion symmetry, the 
topological structure is determined by parity properties at the 
eight high-symmetry points, k*,, in the 3D BZ, which are in- 
variant under time-reversal, up to a reciprocal lattice vector: 
k*, = -k*, H- G. In our case, these symmetries require that 
U//(k) = PU//(-k)P-i and'K,//(k)^ = T9ieff(-k)T-\ 
where the parity matrix P and the unitary part of the time- 
reversal operator T are given by 



P = 



T = 



10-2 



(18) 



where 0-2 is the second Pauli matrix. For any space-inversion- 
odd form factor, it follows immediately that <l>r(k) = at 
a high-symmetry point. Hence, the Hamiltonian at this high 
symmetry point is simply 'Heffik*^) = (^k,-, + ^/k;, )//2 H- (^k;, - 
Efk;„)P/2, where / is the four-dimensional identity matrix. 

The parity at a high symmetry point is thus determined by 
6,„ - sgn(^k;„-£^/k;,)- Four independent Z2 topological indices 
(vo; vi , V2, vs)'^, one strong (vo) and three weak indices (vi_2,3) 
can be constructed from 6m' (i) The strong topological index 

8 

is the product of all eight i5,„'s: Isji - (-1)"° - Yl - ±1; 

m-l 

(ii) by setting kj = (where j - x, y, and z), three high- 
symmetry planes, Pj - |k : kj = 0|, are formed that contain 
four high-symmetry points each. The product of the parities at 



these four points defines the corresponding weak-topological 
index, /" = (-1)''" = H = ±1 (a = 1,2,3) with inte- 

gers corresponding to the axes x, y and z- The existence of the 
three weak topological indices in 3D is related to a Z2 topo- 
logical index for 2D systems (a weak 3D TI is similar to a 
stack of 2D Z2 topological insulators). Because there are three 
independent ways to stack 2D layers to form a 3D system, 
the number of independent weak topological indices is also 
three. A conventional band insulator has all of the four indices 
^STi = ^wTi = ^WTi ^ ^WTi = +1 or equivalently (0;0,0,0). An 
index / = (-1) (v„ = 1) indicates a Z2 topological state with 
the odd number of surface Dirac modes. In a KI the symme- 
try index 6,„ of a particular high symmetry point m is nega- 
tive provided ^k;„ < £/k,*, is lower the f-energy fi/k;,,- Thus if 
^K=o < £/k,*, at the r point, while ^k^^o > Ef^*^ for all other 
symmetry points, then /sti - -1, and hence the Kondo in- 
sulating state is a strong-topological insulator, robust against 
disorde]E2121l Weak-topological insulators and topologically 
trivial insulators can in principle be found for different band 
structures and different values of fi/k;,- A particularly interest- 
ing possibility is to tune topological phase transitions between 
different types of insulators (e.g., by applying a pressure). Al- 
though we have been specifically considering a tight-binding 
model with a primitive unit cell, all our conclusions apply di- 
rectly to systems adiabatically connected to this model. 



I solve the mean-field equations ( 16 1 numerically. For a 
given value of hybridization and temperature I analyze the 



topological structure of the effective Hamiltonian ( 14 1 using 
the prescription outlined above. The results are shown on Fig. 
1 . First I note that when N = 2 the weak topological insulator 
(WTI) state precedes the strong topological insulating state in 
agreement with the earlier studies^'^. For N - A {k - 2) the 
region where WTI exists shrinks and is fully absent for = 6 
{k - 3). It is a quite surprising observation for it is a special 
case when the slave-boson mean-field theory results crucially 
depend on the number of fermionic flavors. In other words, 
here we find an example when there is no adiabatic connec- 
tion between the phases for the realistic case of N - 2 and 
N ^ CO, i.e. when the mean-field theory is exact. The reason 
for the disappearance of WTI phase can be easily traced to the 
condition for the WTI: half of the (5m's must be negative, while 
the other half must be positive. However, as it directly follows 
from the solution of the mean-field equations, this condition 
can never be satisfied when k > 2. Lastly, the results on Fig. 
1 are consistent with the ones which have bee n obt ained for 
the low-energy version of the Anderson modem There it 
was found that WTI is stabilized for /-level energy close to 
the chemical potential for the conductions: this situation cor- 
responds to the /-level occupation ny ^ 1 . With an increase 
of hybridization, n j is reduced as system shows more mixed 
valent behavior, so that nj < I. For N - 2as soon as insulator 
becomes a strong topological insulator, nf ^ 0.8. 

Generally, fluctuations around large-A^ mean-field solu- 
tion introduce interaction between the newly formed heavy- 
fermions. Strictly speaking one needs to prove that fluctua- 
tions of the amplitude and phase of the slave-bosons do not 
break the newly formed state. Due to the presence of the bulk 
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gap, however, we do not expect that fluctuations wifl lead to 
the substantial modification of the ground state. The separate 
issue, of course, is the effects of the fluctuations of the metal- 
lic surface states. Specifically, whether the interactions may 
lead to an opening of gap at the surface. The detailed inves- 
tigation of that problem goes beyond the scope of this paper 
and we leave it for the future studies. 



B. tunneling into topological heavy-fermion semiconductors 

Scanning tunneling spectroscopy (STM) of heavy-fermion 
metals has become an active to pic of experimental and the- 
oretical research in recent yearP^El jjj this regard, an in- 
triguing question is whether the STM can directly probe the 
metallic surface states in a topological heavy-fermion semi- 
conductor In this Section I address this question by evaluating 
the differential tunneling conductance into a weak topologi- 
cal heavy-fermion semiconductor In what follows I will use 
the mean-field theory discussed above for the case of SU(2) 
group, i.e. two flavors of fermions. 

To evaluate the tunneling conductance, we will model a 
bulk system by a stack of L planes along the z-direction and 
diagonalize the Hamiltonian. The resulting Hamiltonian ma- 
trix H„„, has blocks along the diagonal given by (fT4li, which 



has blocks along the diagonal given by 
describe the hopping and hybridization within each plane and 
the off'-diagonal parts describing the hopping and hybridiza- 
tion between the planes. Since the in-plane momentum is a 
good quantum number, the dispersion of the conduction elec- 
trons is given by B{k^,ky) = -(f/4)(cos A;^; + cos ky), while 
the dispersion of the / electrons is described by ef(kx, ky) - 
(f//4)(cos -I- cos A:,,) with tf = O.lf. In addition, we have 
taken the form factor matrix in the form: 



{VisinkjcTx + sinA:, ^",,), within the planes, 
iV^T,, between the planes (upwards), 
-iV^T~, between the planes (downwards). 



(19) 



Within each plane the conduction and /-electrons are de- 
scribed by the four-component sprnor 



/kx,l ''Ik 



x,2 flk^,l ■flk^,2^' 



(20) 



where I labels the layer and = (kx, ky). Lastly, the Hamilto- 
nian describing tunneling between an electrode and a sample 
is 



H,un = T, ^(/7j;.cik^,^ -Hh.c.) 



+ Tf 2 {[<^r(K)]craPlflk,.a + h.C.) 



(21) 



Here we have assumed that the tunneling involves conduction 
and /-electron states on the surface layer only (Z = 1). As we 
will see below, the presence of the form-factor in the second 
term is crucial for the cotunneling events, which ultimately 
give rise to the Fano lineshape for the diff'erential tunneling 
conductance. 




FIG. 2: Plots of the differential tunneling conductance gio)) for the 
stack of L layers of heavy-fermion semiconductors. Panel (a) shows 
g(a>) for the stack of L = 2 layers: for this case there are no states 
in the gap. The band structure as a function of the momentum in 
the 2D BZ is shown on inset. Panel (b) shows g(w) for the stack of 
L = 8 layers corresponding to the weak topological insulator (even 
number of Dirac nodes in the gap). The asymmetry in the tunneling 
conductance is due to the cotunneling processes into conduction and 
/-electron states. 



If we now assume that the tunneling electrode is in equilib- 
rium state with the surface, the tunneling current as a function 
of the voltage, /(V), is 



KV) 



2e 



oo 



d(jjp,ip{(jj — eV)[nf(a> — eV) — n/r (w)] 



^ J ■ (22) 

— OO 

X Im [\TfG,,(co) + \T/Gff(oj) + 2\TA\Tf\G,f(oj)] , 

where p,,,, is the STM tip density of states (DOS), n/r(w) is the 
Fermi distribution function and Gai,{a)) are the advanced local 
single particle Green functions: 



^ , , V C,/ki)^'^.^(k^) 

Gai(w) = > -r— — . 

w - e^(kj_) - 10 



(23) 



A,k± 



In the expression above, £^(kj^), 0a„i(kj_) denotes the set of A 
eigenvalues and the corresponding eigenfunctions I obtained 
by diagonalizing the Hamiltonian //,„„, while subscripts a, b 



refer to the components of the spinor ( 20 1 at the surface (Z = 
1). In real materials, however, electronic correlations as well 
as disorder lead the broadening of the /-electron level. One 
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way to account for these effects in the differential tunneling 
conductance g(V) = dl/dV is to consider the complex valued 
quasiparticle energies: 



with the quasiparticle width given hy ^'^ 



< T* 



|£,kj/[l + log(|£,kJ/r^)]^ £.!k, > T*. 



(24) 



(25) 



Here T* is the temperature corresponding to the opening of 
the hybridization gap, or the temperature at which the first 
non-trivial solution of the mean-field equations appears. In 
SmBe, for example, T* ^ lOOK. With these provisions, we 
evaluate the differential tunneling conductance for the set of 
the parameters corresponding to the weak-topological insula- 
tor in 3D translationally invariant system. I show the results 
on Fig. 2 for a fully gapped states (top panel) and metallic 
state (bottom panel) corresponding to the weak topological 
insulator. Both curves have characteristic asymmetry due to 
the co-tunneling processes into conduction and /-states. The 
finite value of the g{V) are zero bias are due to the finite width 
of quasiparticle states, Eq. (25 1. From our results we see 
that gwTiiV ~ 0) in the case of weak topological insulator is 
significantly enhanced in comparison with gBi(V ~ 0) (band 
insulator), which is not surprising. In addition, gwniV) shows 
higher asymmetry then gsiiV). Nevertheless, it is seems to be 
a quite challenging task to argue in favor of the topologically 
protected surface states solely on the STM data. 



III. CONCLUSIONS AND OUTLOOK 



i.e. cubic, point group symmetry. In this case, the analysis 
of the crystalline field split /-ion multiplets for the valence 
configurations corresponding to the total angular momentum 
J - 5/2 or J = 7/2 shows that only N - A degenerate mul- 
tiplets can give rise to the nodeless hybridization gap. Such 
a scenario can be realized in heavy-fermion semiconductor 
SmBg. Indeed, finite metallic conductivity below T ^ 5K 
may serve as a signature for topologically protected metallic 
surface states. The fact that conductivity grows with the im- 
provement of the quality of the samples qualitatively supports 
this idea. Indeed, recent theoretical works^^ have explicitly 
demonstrated that the presence of relatively strong disorder on 
the surface of a strong topological insulator will substantially 
disrupt these states. The physical reason for such a behavior 
is that the impurity induced states propagate well below lead- 
ing to the diffusive behavior in the surface transport. In that 
regard, when the strength of disorder potential is comparable 
to the bulk gap, the 2D Dirac theory description of the 3D 
topological insulators is not valid''^. 

An important issue for the subsequent study is a role of 
fluctuations around the mean-field solution. For the band in- 
sulator, one may argue that the fluctuations effects (generally 
of the order of 1/A^) do not lead to any significant changes 
providing only small corrections to the gap itself. Situation 
becomes drastically different when the metallic surface states 
are present. In particular, fluctuations in slave-boson fields 
lead to the effective interactions between the conduction and 
/-electrons, which in principle may lead to the opening of the 
gap at the surface as well. The detailed investigation of these 
effects goes beyond the scope of this paper and I leave it for 
the future. 



In this paper I have analyzed the low-temperature properties 
of the generic heavy-fermion semiconductors using the large- 

slave-boson theory. Specifically, I have provided an evi- 
dence for the formation of the topologically non-trivial elec- 
tronic states at the surface of these materials. The results re- 
ported in this paper are in agreement with those obtained us- 
ing different approach based on the low-energy analysis of the 
Anderson lattice model. The phase diagram obtained within 
the mean-field analysis implies that the strong topological in- 
sulating phase is likely to be observed in materials with high. 
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